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Abstract 

We describe a gradient flow structure for the inviscid Langmuir layer 
Stokesian subfluid model introduced recently by Alexander et al [1]. 

1 Introduction 

In a recent article Alexander et al [1] introduced a model for the evolution of 
molecularly thin Langmuir layers on the surface of a subfluid. The Langmuir 
layer consists of domains of different phases, and the simplest model is a study 
of layers with only two phases. For a quiescent subfluid, the evolution of the 
Langmuir layer is driven by a line tension at the interface between phases and 
damped by a Stokes flow in the subfluid. The resulting free boundary problem 
is called the inviscid Langmuir layer Stokesian subfluid (ILLSS) model. 

In this brief note we comment on the mathematical structure of the ILLSS 
model. Alexander et al showed that the perimeter of the interface determines a 
natural energy of the system which is dissipated by the Stokesian subfluid. We 
show that this energy-energy dissipation relation arises from a gradient flow. 

This observation is based on a comparison between the ILLSS model and 
the Hele-Shaw flow. Alexander et al showed that the ILLSS model admits a 
boundary integral formulation similar to that of the Hele-Shaw flow, and it 
is known that the Hele-Shaw flow has a gradient structure. The first such 
formulation is attributed to Fife in [2], We follow a distinct approach due to 
Otto [4]. The ingredients of a gradient flow are a manifold (the phase space) 
with a metric and an energy functional. We find that as in the Hele-Shaw 
flow, the energy in the ILLSS model is the perimeter of domains. However, it 
differs from the Hele-Shaw flow in the choice of manifold and metric. The main 
heuristic idea is that the metric must reflect the dissipation, and here it occurs 
in the Stokes flow in the subfluid. 
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There is also a striking difference with the Hele-Shaw flow. An interest- 
ing feature of the ILLSS model is the remarkable stability of slender domains 
(tethers), both in computations and experiments. This is in sharp contrast 
with computations that suggest singularity formation for tethers evolving by 
the Hele-Shaw flow [2]. Our observation may be of value in addressing these 
harder questions. 

2 The model 

The ILLSS model is as follows [1, §2.4]. The reader is referred to [1] for a 
description of the assumptions that underlie the model. 

The subfluid occupies the semi- infinite domain B — {(x,y, z) G IR 3 |z < 0}. 
The velocity u = wi + uj+wk and pressure P satisfy the Stokes equations 

V-u = 0, z<0, (1) 
Au = VP, z < 0. (2) 

The normal velocity w vanishes on the boundary z = 0. We assume that u and 
its derivatives decay sufficiently rapidly as x 2 + y 2 + z 2 — >• oo. 

The Langmuir layer is modeled as the surface z = of B. It decomposes 
into two Langmuir layer domains f2 and fl c separated by a moving boundary 
<9f2. The motion of the boundary is determined by force balance and kinematic 
conditions. 

The force balance is as follows. The Langmuir layer z — is assumed to be 
in hydrostatic equilibrium. Then the surface pressure LI on z = 0, the tangential 
surface stress F s , and the line tension F; satisfy 

V±II = F a +F ; , (3) 

where Vj_ = d x i + d y j. The Stokes flow in the subfluid determines 

F s (x,y) = -u z (x,y,0)i-v z (x,y,0)j. (4) 

The surface pressure LI jumps across the domain boundary <9f2 because of the 
curvature and the line tension. This may be written 

F ; = KnJ(d) (5) 

where k is the curvature of 9f2, n is the unit outward normal, S is the Dirac 
delta, and d denotes the normal distance from <9f2. The convention in [1] is that 
k < for convex J7. 

We now consider the kinematic conditions. The surface velocity of the Lang- 
muir layer is denoted U. We assume the Langmuir layer is incompressible 

V±-U = 0, (6) 

and that the surface velocity is compatible with the subfluid velocity 

V(x,y)=u{x,y,0). (7) 
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Finally, a point T on the boundary is advected by the surface velocity 



^=U. (8) 

Dt v ' 

This completes the specification of the ILLSS model. 

Alexander et al assume in addition that the flow field leaves horizontal sec- 
tions invariant. As a consequence the vertical velocity w = vanishes in B (not 
just on z = 0). This assumption was made in earlier work [3, 6] and allows 
considerable simplification via the introduction of streamfunctions. 



3 Gradient flow structure 
3.1 The framework 

The ingredients of a gradient flow structure are a manifold M. with a metric g 
and an energy functional £. If m e M, the gradient flow may be written in 
weak form as 

g(rh,v) = -(d£(m),v), veT m M. (9) 

Here (•, •) denote the duality pairing between 1-forms (elements of T m A4*) and 
vectors (elements of T m A4). The role of the metric is to convert the 1-form 
—d£{m) into a vector rh € T m A4. 

We now show that the ILLSS model is a gradient flow. The calculations 
are similar in spirit to [4]. The manifold Ai consists of C°° diffcomorphisms 
of B that preserve orientation, volume and the surface area on the boundary 
dB. We denote a typical element of M. by <p and its restriction to dB by $. 
Diffcomorphisms in M. correspond to flows generated by divergence free vector 
fields in B with a divergence free trace on dB. Assume v is a C°° vector field 
B — >• M 3 with trace V(x,y) = v(x,y,0) and V • k = 0. The associated flow 
<P v (t) is the solution to 

a r ^ v (T) = vo^ v (r), p(0)=ld. (10) 

This also implies 

a T $ v (r) = Vo$ v (r), $(0)=ld. (11) 

These calculations allow us to identify the tangent space T V M. We first define 
the linear space of Eulerian velocity fields that form Ti^M. 

V={vgC oc (B,M 3 )| V-v = 0, V_l-V = 0, V-k = 0}. (12) 

Let v,p denote the composition v o ip. Then the tangent space at ip € M. is 

T v M = {v v :=vo<p,v £V}. (13) 

We now introduce the metric. Let ip G M and v e V. We define 

3 

SVK» v v)= / Tr (Vv T Vv) dx = V / {d^vtf dx. (14) 
Jb iJ=1 Jb 
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If the integral vanishes, v is a constant that vanishes because of the boundary 
condition. Thus, g is indeed a metric. Note also that it is enough to specify the 
metric by (14). This is equivalent to 

3 v (v v ,w v )= / Tr (Vv T Vw) dx, v,weV. (15) 

J B 

The principal heuristic idea is that the metric is determined by dissipation. All 
the dissipation in the ILLSS model is in the subfluid. The expression (14) is 
precisely the dissipation in a Stokes flow in the subfluid. 



3.2 The unconstrained ILLSS model as a gradient flow 

The natural energy of the ILLSS model is the perimeter of the interface <9S! t . 
If an initial phase configuration ilo C dB is fixed, every ip e M defines a new 
configuration $(f2o) (recall that $ is the restriction of ip to dB). We set 

£(<p) = Per($(n )). (16) 

We are now ready to formulate the ILLSS model as a gradient flow. Equation 
(9) takes the form 

g v (dtip,v v ) = -{d£{y),v v ), v^eT^M, (17) 

where the flow of diffeomorphisms ip(t) is generated by the 'true' velocity 

d t tp = uo<p, a ( $ = Uo$, <p(0) = Id, $(0)=ld. (18) 

We show that equations (18) and (17) are equivalent to the free boundary evo- 
lution of Section 2. 

We first compute the left hand side of (17). We use equation (15), integrate 
by parts and use (4) to obtain 

g v (dt<P,v v )= [ Tr (Vu T Vv) 

J B 

= - Au(x) ■ v(x) dx+ u z V 
J b Job 

= -/ Au(x) • v(x) dx - I F s -v. (19) 
J b Job 

The right hand side of (17) is computed as follows. The first variation of the 
perimeter is given by the curvature [5, §2.9]. That is, 

(^),v)=^^ 



dr 



T = 



[ nn ■ V = - f nS(d)n ■ V = - / F; • V. (20) 

JdQ JdB JdB 
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(The —sign follows from the convention of [1] that the curvature of a convex 
domain is negative). We combine (19) and (20) to see that (17) is equivalent to 

- [ Au(x) ■ v{x) dx = [ (F ; +F s )-V, vgV, (21) 
Jb Job 

with Fi and F s defined by the constitutive relations (4) and (5). It follows that 

- f Au(i) • v(x) dx = 0, (22) 

JB 

for every v G V with V = 0. 

The Helmholtz decomposition for a vector field f G L 2 (B) is 

f = VP + w, V-w = 0, W-k = 0. (23) 

The pressure P is obtained as a solution to the Neumann problem 

AP = 0, z < 0, d z P = f • k, z = 0. (24) 

P is arbitrary up to a constant, but VP G L 2 (B) is unique, as is w. We apply 
this decomposition to f = Aw. Then (22) yields 

- / w(x) • v(ar) dx = 0, (25) 

JB 

for every v G V with V = 0. It follows that w = 0. Thus, u solves Stokes 
equations (l)-(2). Further, since Au = VP, we integrate by parts in (21), and 
use V • k = to find 

f (F i + F s )-V = 0, (26) 

JdB 

for every smooth, divergence free vector field V on dB. We apply the Helmholtz 
decomposition again (but now in IR 2 ) to deduce the existence of a surface pres- 
sure II such that (3) holds. Finally, the kinematic condition (8) is immediately 
implied by the fact that solves (11) with V = U. This establishes the 
equivalence between (17)-(18) and the ILLSS model of Section 2. 

3.3 The constrained ILLSS model as a gradient flow 

We now consider the constrained flow considered in detail in [1]. The constraint 
is that the flow leaves horizontal planes z ^constant invariant. The flow is 
now restricted to the submanifold Af C M. that consists of ip € A4 such that 
ip(x,y,z) = (ifi(x,y, z),if2(x,y, z), z). The corresponding tangent space T V N 
consists of vector fields v o ip as in (13) with the additional constraint v • k = 0. 
We now claim that the constrained ILLSS flow of [1] may be written as the 
gradient flow 

g v {d t <p,v v ) - -(d£(ip),v v ), v v G T V M, (27) 

where ip(t) also solves (18). Again, we must show that u satisfies the equations 
of Section 2. 
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The left and right hand sides of (17) are computed exactly as earlier and yield 
(21) as earlier along with the constitutive relations (4) and (5). We again choose 
v with V = to see that (22) holds for every v that is divergence free and satifies 
v • k = and V = 0. The argument involving the Helmholtz decomposition of 
Au is modified as follows. We now decompose a vector field f e L 2 (B) with 
f • k = in B as in (23) with the additional constraint that w • k = in B. 
It now follows that P = 0, since the Neumann boundary condition in (24) is 
simply u • k = 0. We apply this decomposition with f = Au to (22) and obtain 
(26) and the conclusion w = as above. This shows that the velocity field u 
is harmonic (thus also a solution of the Stokes equation with P = 0). It now 
follows immediately from (21) that (26) holds for every divergence free vector 
field V. This yields (3). This shows the equivalence between (27) and (18) and 
the constrained model of [1]. The conclusion that u is harmonic is deduced 
in [1] using streamfunctions. 
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